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differential equations is presented. Consequently, the real number set gets replaced by an 
extended physical set, each element of which is endowed with an equivalence class of 
infinitesimally separated neighbours in the form of random fluctuations. We show how 
a sense of time and evolution is intrinsically defined by the infinite continued fraction of 
the golden mean irrational number g"" 1 , which plays a key role in this extended SL(2,R) 
formalism of Calculus. Time may thereby undergo random inversions generating well 
defined random scales, thus allowing a dynamical system to evolve self similarly over 
the set of multiple scales. The late time stochastic fluctuations of a dynamical system 
enjoys the generic 1/f spectrum. A universal form of the related probability density is 
also derived. We prove that the golden mean number is intrinsically random, letting all 
measurements in the physical universe fundamentally uncertain. The present analysis 
offers an explanation of the universal occurrence of the golden mean in diverse natural 
and biological processes. 
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1 Introduction 



In Calculus, a variable changes by ordinary shift operations (translations). Every point in the 
real axis is conceived as a structureless point. In a dynamical problem, the time evolution 
of a system occurs over a finite number of characteristic scales as inherited from the under- 
lying differential equation. The generation of fractal- like self similar structures over multiple 
(dynamically generated) scales, either in the phase space or purely in the time sector, needs 
explicit nonlinearity at the level of the governing equation. In the following, we present a 
new class of scale free stochastic solutions of the simplest ordinary differential equation. A 
representative solution of the new class, which is exact and finitely differentiable, defines an 
extension of the real number set, endowing every number with a nontrivial neighbourhood of 
fluctuations. These fluctuations experience an intrinsic universal evolution, which, in turn, 
generates a generic irreversible sense of time defined by the infinite continued fraction repre- 
sentation of the 'golden mean' irrational number. Consequently, time may undergo random 
inversions, allowing a dynamical system to explore stochastically, new scales, usually unavail- 
able in the ordinary dynamics. The corresponding evolution of the system would, therefore, 
enjoy a generic late time stochastic fluctuations with 1/ / spectrum [1]. We also obtain an exact 
form of the universal probability density, recently found [2] to occur naturally in a wide class 
of scale free nonlinear processes in Nature. These exact results provide a natural explanation 
of the universal occurrence of the golden mean in diverse natural and biological processes. 

The relevance of the present work may also be seen in the light of the recent works [3,4] 
uncovering new relationships between time and the number theory. In ref [3], the evidence 
of 1/f noise and scale free self similarity in the prime number distribution is pointed out. 
In ref [4], Planat reported the observation of discrete time jumps in beat frequencies in the 
context of a superheterodyne receiver. These random jumps are shown to generate 1/f noise 
in the oscillator frequencies. Further, the origin of 1/f spectrum is related to the arithmetical 
summatory functions such as the Mobius and the Mangoldt functions which arise in connection 
with the Reimann zeta function. Sometime back, Robinson, on the other-hand, developed [5] a 
theory of infinitesimals on the basis of mathematical logic, thereby giving rise to a valid model 
of a nonstandard extension of the standard framework of analysis. As it turns out, the present 
extension of Calculus is a realization of the nonstandard analysis, in which infinitesimally small 
numbers are shown to have dynamical properties. We note also that our exact results seem to 
provide a correct mathematical framework for the theories of fractal space-time being developed 
by several authors [6,7,8] . A significant part of El Naschie's work [8], in particular, explores the 
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significance of the golden mean in Cantorian fractal spacetime and other related branches of 
Physics and Mathematics. A detailed comparison of the present results with those of Refs. [6-8] 
would require separate investigations. Incidentally, the possibility of fractal solutions to linear 
differential equations seems to have been pointed out first by Nottale [6].[] 

To indicate how the exact results, reported here, are obtained, starting from a heuristic 
definition of time inversion, we present the exact results in Sec. 4. In the preceding two sections 
we indicate how the exact class of solutions is derived in Sec. 3, using an approximate analysis, 
based on a local definition of time inversion, explained in Sec. 2. Applications of the results 
are discussed in Sec. 5. We close the presentation with some further remarks pointing out the 
future scope of this extended formalism of Calculus. 

2 Time inversion 

Let us consider the simplest linear dynamical system given by 

dJ = (! + Kt ) x C 1 ) 

where t denotes the dimensionless time ( we scale t suitably to adjust the dominant scale 

of evolution to t ~ 1) and k is a small, slowly varying (almost constant), parameter. In 

the ordinary Calculus, this equation has the 'standard' solution x s oc exp(t + |k£ 2 ), with 

no (self-similar) fluctuations, unless the equation is explicitly nonlinear, for instance, through 

K = n(t, x, x, x, . . .). Our intention is to show that eq (HD can accommodate a new class of 

nonlinear stochastic solutions, even for a constant k, under a simple, but general assumption 

that time may change by inversions as well. 

Let us recall that ordinarily a change in time, in the vicinity of a given instant t is indicated 

by a pure translation t = t + 1 = to + (t — to)- 111 fact, the last equality is an identity (valid 

for all t). By an inversion, on the other hand, we mean the following. Let t± denote times 

t J$ 1 and t ~ 1 respectively (from now on t denotes the rescaled variable t — > ~). Then 

close to t — 1, the inversion t- = 1/(1 + (£+ — 1)) leads to the constraint 1 — £_ = t + — 1. 

The parametric representation of inversely related times is obviously given by £_ = 1 — t and 

t + = 1 + i, < i << 1 (so that the constraint reduces to an identity, valid close to t — 1). 

With this reinterpretation, time inversion in the vicinity of t = 1 assumes a form analogous to 

a pure translation. If translation is considered to be the most natural mode of time increment, 

then there is no compelling reason of ignoring inversion as yet another natural mode of doing 

1 Prof. El Naschie has pointed out that it is in-fact G. N. Ord who first conceived the possibility of fractal 
solutions to linear differential equations. 
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this. Consequently, it seems reasonable to assume that time may change from t_ to t + not only 
by ordinary translation over the period t + — i_ = 2i, but also instantaneously by an inversion. 
Let us clarify the physical content of the definition further through several remarks. 

1. Note that the above definition gives a new nontrivial solution t- = l/t+ to the constraint 
t_ +t + = 2, in the vicinity of i = 1 over the linear solution t- = 2 — 1 + , ordinarily thought to be 
the only possible solution. In view of this nonlinear possibility, the change (flow) of time could 
be visualised as an SL(2,R) group action, when the nontrivial SL(2,R) action is realized only 
in a neighbourhood of a point, t — 1, say. One may thus imagine that time flows, for example, 
from t = to t = t- by translation, and then may switch over to t+ by inversion t + = for 
another period of linear flow etc. In the next sections, we show how this second and subsequent 
periods of linear flows are actually realized over scales of the form t n = u n t, n = 1, 2, . . ., where 
v = = 0.618033939..., is the golden mean irrational number.^ 

2. The definition of time inversion has an inbuilt uncertainty. The moments t± are not 
well defined, except for the fact that these be close to t — 1, thus elevating i, and hence time 
itself, to the status of a random variable. Note that this sort of uncertainty (randomness) is not 
allowed in the framework of linear Newtonian time. In the present SL(2,R) framework, time 
may, therefore, undergo small scale random fluctuations in the neighbourhood of every instant 
(t << 1). Two instants t± joined by an inversion, can therefore be in a continual process 
of transport between each other by randomly flipping their signs. The possibility of random 
flipping also indicates that there is an inherent uncertainty in the actual determination of an 
instant close to t = 1. In case one claims that he/she is in the moment i_, then there is a 50% 
chance that he/she is actually at the moment t + . From now on, we denote by t this stochastic 
behaviour of time, when the ordinary Newtonian time is denoted by 77. 

3. The ordinary nonrandom variable (time) 77, that we are accustomed to at moderate scales 

should be retrievable from the stochastic t in the mean 77 =< t > (upto a rescaling). Physically, 

it means that the small scale fluctuations near every point of t cancel each other in such a 

way to yield the average coarse-grained time sense 77 to our experience. These fluctuations 

would, however, become important to determine the small scale structure of time and hence 

of a dynamical system. By inversion, these would also have nontrivial influences on the long 

time behaviour of the system. Incidentally, we note that the formal definition of inversion is 

still valid for 77. For a nonzero mean i, for instance, inversion between r)± — 1± < i > could 

2 The standard notation of the golden mean in the literature is <fi. Our symbol v in this work is motivated 
by the fact that the golden mean appears here as a universal scaling factor of time. The symbol <f> is more 
often used to denote a quantum state function as well. In Sec. 4-5, we, however, denote the stochastic time-like 
character of the golden mean by cf>. 
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then be implemented at a well determined instant. For a zero mean i(< t >= 0) , on the 
other hand, the inversion constraint collapses to r] = 1. Thus allowing for a random, zero mean 
i is equivalent to probing small scale structure of a point, t — 1, where by point we mean a 
'physical point' which is meaningful only in the context of an accuracy limit (level of resolution), 
in contrast to an ideal structureless point . As an example, let us consider a situation when the 
measurement of a duration (an interval) is allowed only upto the first decimal accuracy. In that 
case, points separated by a distance less than .05 are indistinguishable, and should be treated 
as an equivalence class. In this sense, the point t — 1 is an equivalent class, where the symbol 
'1' is only a convenient representative of the class. In the ideal case of infinite accuracy, the 
class '1' has only one element in the ordinary (linear) Calculus. In the presence of inversions, 
it now follows that even an ideal point ( in the limit of infinite accuracy) has infinite number of 
nontrivial members separated by infinitesimally small, zero mean random scales. Since i(t) is 
a positive, monotonically increasing stochastic variable, the nontrivial equivalence class of an 
ideal point has, at the least, the cardinality of continuum. In fact, it can be higher, as is shown 
in Sec. 4. (It will be evident that zero mean is not a necessary restriction.) 

4. It turns out that the concepts of inversion and stochastic time are deeply related to the 
(practical ) limitation on the exact measurability of a duration. Let, for instance, the moments 
t_, t + belong to the equivalence class of the 'physical point' t — 1. Then for any practical 
purpose, t- — t + — 1. That means, in turn, that any measurement over a finite period of time 
would fail to make any distinction between t_, t + and 1. The associated fluctuations between 
t- , t + which remain imperceptible over a pretty long time can, however, act as potential seeds 
for self-similar evolutions over scales t n = u n t, which can grow slowly to interfere with each 
other, thereby deflecting the standard nonrandom evolution, to a universal pattern of stochastic 
fluctuations (c.f.,Sec.3). 

5. It is worth comparing the present definition of a local time inversion with the usual 
(global) time reversal (inversion) symmetry of an equation of the form eq(l). The usual time 
reversal symmetry means that the system x(t) evolves not only forward in time from t to 
t + h, h > 0, but it can also evolve backward; i.e., the state x(t) can be reconstructed from the 
state x(t + h). The parameter t in eq(l) is thus 'non-directed', giving rise to the problem of time 
asymmetry. As remarked already, the new class of solutions indicates an inherent irreversibility 
in the time sense. 
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3 Fractal Solution 



In the framework of SL(2,R) stochastic time, eq (|l|), rewritten as 



dx — (1 + nt)x dt 



(2) 



assumes the status of a stochastic differential equation. We now present a treatment of the 
above equation when the variable t is assumed to act as the ordinary 'coarse grained' time r], 
except in the vicinity of t — 1 (the dominant scale of evolution). Thus over the period (0,1), 
the system evolves along the standard evolutionary path x = ~ e v (k being small). The 
small scale fluctuations in the system variable x, as inherited from the fluctuations in t would 
remain un-noticeable in comparison with this mean evolution. One can, however, probe these 
fluctuations, provided the mean evolution xq is removed from the actual (total) system x via 
the ansatz x = xox±, where x\ is the purely 'fluctuating' component satisfying the reduced 
equation 



To prove the self similarity of the small scale fluctuations, and to see how these grow and 
influence the late time properties of the mean evolution, we probe the neighbourhood of t = 1 
more and more closely, using scales t n , separating in each step the relevant residual evolution 
into 'mean' and 'fluctuation' respectively. To see this explicitly, let us assume, for definiteness, 
that the actual removal of the mean in eq(||]) is accomplished at an instant t- ^ 1. We note 
, however, that the specification of the exact moment is physically impossible. Thus the moment 
t- must have an inherent uncertainty, which we represent by an ansatz t- — 1 — rfj, where, 
r is a discrete variate with moments < r n >^ Tq, < r >= ro, so that < £_ >=?]_. For the 
sake of clarity, one may also assume that t + and t- belong to the physical equivalence class of 
l(c.f., remark 4), when a point is defined correct to m decimals (so that the accuracy limit is 
10 _m ). Because of the random inversions t + <-» t_ in the time coordinate, the residual system 
pair Xi(t^) and x\{t + ) also undergo fluctuations between them. However, these fluctuations 
would remain unobservable at least over a sufficiently long period of the form < i] < 1 + a, 
provided ro is sufficiently small and a ~ l/(10"Yo). Now utilising the golden mean partition of 
unity: z/ 2 + z/ = l, z/>0, one can realise a scale changing SL(2,R) transformation 



dx 



i — 



Kt X\dt 



(3) 



t + = 1 + rrj 



1 + fv r\\ 



(4) 



1 — rrji 
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where r]x = r ui] = v{r] + — 1) << 1, r = r/r . Using eq(^) and the inversion constraint 
dt_ = —dt + = — rdf], eq(|3|) gets transformed to 

— dxi = fA(l + rvT\) xidTi (5) 

where T\ = In 771, rji = 1 + 771 = 1 + ^(77+ — 1) and X = k/u, when we make use of ln(l + 7771) « 
rfji ~ rT\. This equation, valid close to fj — 0, (77 = 1), and self similar to eq(0) , describes 
the small scale evolution of the first generation fluctuation x±. (One can, indeed, recast eq(^) 
exactly to the form eq(^) m the stochastic time variable t\ = fTi.) The (-) sign is a signature of 
inversion. Note that this also avoids the possibility of a backward flow of time at the expense of 
deflecting the direction of system evolution. As time 77 continues to flow from 77+ ~ 1 onwards, 
dragging fji along with, the small scale fluctuation X\ also gets amplified following eq(|5]) and 
assumes the status of the original system variable x over time 77 ~ 0(e(l + v) — i>), when a 
second generation transition to the scale T 2 becomes permissible. Note that the self similarity 
of eq([5]), relative to the time variable t\, with eq(|l|) tells that T\ would act as the ordinary time 
for Ti e (0, 1) for the evolution of X\. Factoring x\ = e~ xtl X2, one thus gets the 2nd generation 
replica of eq(|l]) for x 2 i n the time variable T 2 and hence this method of self replication over 
scales T n could continue ad infinitum. An infinite string of iterations, as above, thus leads to a 
new solution of eq(|I]) in the form 

x oc e v-FW-V>+...) (6 ) 

Note that T 2 = lnr/ 2 , r] 2 = l+v(Ti— 1) = l + z/(ln(l+z/(r/ — 1)) — 1)) and hence the nth generation 
scale T n is related to rj by n nested natural logarithms. Consequently, the over all fluctuations in 
the system Xf oc x/xq = ( ^^ ,,, ) rA , incorporating influences of all scales, has the asymptotic form 
x f - ( ln ;!nln^- ) r " A = ^" r> , 3 as rj - 00. Here, the exponent /x = A(l - gj), a = is 
a slowly varying function of rj. A few remarks are in order. 

1. The choice of the moment t = 1, around which the evolution is probed, is for the sake of 
convenience. For any point t G (0, 1), the analysis proceeds with the rescaling t — > t/t , k — > 
t^n in eq(|^). 

2. Because of the nested logarithms, contributions from higher order scales are felt slower 
and slower. As stated above the first sign of fluctuation is surfaced only if the system is allowed 
to evolve over a period 77/ ~ 1 + 10 _m r _1 . As an example, in a 'universe' where only the first 
order accuracy (m = 1) is allowed, the fluctuation is first noticed around rjf ~ 2 for a tq ~ 0.1. 
More generally, for a tq = l0~ ( - m+s \ r/f could be arbitrarily large, for a large value of s, even in 
the limit m —>■ 00 of infinite accuracy. 
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3. It is natural to interprete the new solution eq(|6|), with self similar fluctuations over the 
mean solution (random) fractal solution of the equation 



under time inversions. This is, however, in contrast to our aim which was to find such a solution 
for eq([3]), so that stochastic fluctuations would have been around the standard solution x s . In 
the following, we resolve this dichotomy in a more general framework of fractal time, where the 
parameter A is identified with an 'apriori' scale factor of time, rather than a system variable. 
We also show how the correct fluctuation pattern is obtained when A acts as a true system 
variable. 



4 Fractal time 

We have shown how a stochastic behaviour could be injected to time via inversion, while time 
inversion is interpreted as a consequence of the practical limitation of exact measurability 
of a duration. ( Stated more precisely, the possibility of time inversion raises the practical 
measurement limitation to the level of a theoretical principle.) The stochastic nature of time 
already endows time with fractal -like characteristics. To investigate fractal properties in more 
details, let us begin by writing an ansatz for the 'physical' fractal time t as an implicit random 
function of the ordinary time r\ : t±{rj) = f](l ± K'qt^'q^ 1 )) . Here, K > stands for an arbitrary, 
small but random, parameter. The nontrivial second factor would be responsible in determining 
the small scale structure of the physical time in the neighbourhood of every point rj. The choice 
of sign =p in t in the factor indicates explicitly the possibility of inversion near rj = 1 (see below). 
Note also that t± mimics the notation of Sec. 2-3, thereby splitting every point of rj axis into 
an equivalence class {t±} of infinite number of finely separated points. Since each member of 
the class is a function of rj, it has, at least, the cardinality of continuum (c.f., remark 3, Sec. 2). 
We show below that the actual cardinality is 2 C . 

As it turns out, the ansatz represents a new class of exact, stochastic solutions to the 
equation ^ = 1. To verify this, we note (suppressing the distinction temporarily) that, by 
symmetry, both t/rj and rjt, t = t(^ _1 ) satisfy coupled equations of the form a = 1 + k/3 
and (3 = 1 + na, hence t(r))/r) = r}t(r)) for all t. Noting that ^ = ~ ?7~ 2 ^pT, we get || = 
(1 + Ki]t ) + rjKt - K^pr, so that - |) + ft(^pr — rfi ) = 0. It thus follows, k > being 
arbitrary, that 
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d?7 

which is nothing but the desired equation in logarithmic variables. Note that for a nonrandom 
real parameter k, one retrieves the standard solution t — (1 — k)^ 1 ^. For a random k, which 
arises naturally in the context of the present 'local' definition of inversion, we, however, get a 
new random (fractal) solution, which matches (approximately) with the standard linear solution 
only in the mean. We emphasise that the new solution is an exact solution of eq(|8|). Because 
of its inherent scale-free nature, the solution must possess nontrivial fractal characteristics. 
Indeed, eq(||) tells that hit, and hence t/i], must be a function of In 77, so that t = rjcpilnr]). 
Here, 0(ln77) = c(l + K0(ln77 _1 )), and represents a nontrivial (random fractal ) solution of 
^ = 0, c being a real constant. 

Let us note that a straightforward iteration of the ansatz in the form t/i] = l + K + K 2 + ... 
would be, in general, misleading because this geometric series in k apparently hides the slow 
time dependence that is always present in any finite approximants of this infinite series. We 
have already shown how such a slow, residual time dependence in the nth approximant S n = 
1 + k + k 2 + . . . + K, n r]t can influence the dynamics over time 77 ~< k > _n , because of local time 
inversions. 

To explore the role of golden mean u, ab-initio in the present context, let us now reintroduce 
signs '±' to distinguish the variables t±. Let t + /r] = 1 + rr/it^^ 1 ) and t_/r] = 1 — rr)it + (r)i), 
where 7/1 = kr], k = rk, k is an ordinary constant, and r (~ 1) is a random variable, analogous 
to one in Sec. 3. Note that the present form is slightly general from above, but, nevertheless, 
solves eq([8]), because of its scale free nature. This scale free property now tells that the limit 
of t + (r])/r] = 0(77) as 77 — > r/ , < r] < 00 is independent of r] . As a consequence, 4>(rj) is 
a universal random function, defined in the vicinity of 77 = 1, so that t + = r]<f>(rj), 4>(r)) = 
1 + r0(r7jf 1 ). Note that for a sufficiently small k, the variation of 0, which always remains of 
the order 0(1), is very small. Note that £+1(77) = kt + (ri) = t + (rji), by definition. 

Let us now recall that 771 = 1 is an 'ideal' point in the ordinary time axis, in contrast to 
the 'physical point' t — 1, an equivalence class of members of the form {t±} in the geometrical 
axis of the physical time. Now, as 77 approaches 00 crossing A; -1 , the rescaled variable 771 
crosses 771 = 1, running over points such as 77i_ = 1 — a to — 1 + a. An ordinary time 
inversion 77i_ = 77^ (c.f., remark 3, Sec. 2), now induces a random inversion in the physical time 
t-(r)i_) — ► £-(771+) : rt_(?7i + ) = 77^ + /t + (?7i + ), so that r<j>{j]^ ) = l/(f>(r)i). Note that the exact 
moment of inversion is uncertain because of the inherent randomness in the physical time due 
to the r.v. r. Consequently, one obtains t + /rj — 1 + 771 /t + (771), with rji = 1 + a, a — > 0. It thus 
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follows that t + = (jyq, 4>{r]) = 1 + v, which is true not only for rj — > oo, but, in fact, for any r\ 
because of the universality of </>. 

To continue further, let us note that the above result, although proves the unique role of the 
golden mean number in the framework of fractal time, also presents us with a riddle. Apparently, 
one would like to conclude that this solution reproduces the standard solution t = cq of eq(]8|). 
However, the emergence of this unique special value v is unclear in the ordinary framework. 
Recall that the variables t± are intrinsically random. Further, no where in the above analysis 
we have taken mean values to eliminate the underlying randomness. The only reasonable 
conclusion would therefore be the following: the golden mean number <ft represents a universal 
random fractal function which is responsible for small scale random fluctuations in the physical 
time t.0 

To explain the assertion in detail, we need to proceed in a number of steps, establishing a 
number of key results. 

Let us begin by noting that the ordinary solution of eq([8[) defines a 1-1 ( the identity) 
mapping R — > R of the ordinary real set. The new class of random scale free solutions now 
defines an extension of the ordinary real set to the 'physical ' real set P, say. The extended 
solution space of linear ordinary differential equations now consists of solutions which are (i) 
infinitely different iable and (ii) those which are finitely differentiable, stochastic functions. We 
shall verify shortly that the new implicitly defined functions are first order differentiable, with 
discontinuous second derivatives, at the moments when rj changes by inversions. Because of this 
extension, it is natural to expect that the physical set P contains new members not available to 
R. To show that P indeed has nontrivial numbers, let us first distinguish a physical number K p 
from an ordinary positive real number K, where K p = K + N, N being the neighbourhood of 1 
consisting of random physical numbers (fluctuations). Thus, a physical number K p , representing 
a nontrivial equivalence class, is a non-singleton subset of R+, the set of nonzero positive reals. 
Further, by definition K p = K<f>{rj). 

Let P(R) denote the power set of R = R + U (—1,0), P, being the corresponding extension, 

and N the set of discrete subsets in R. Then the cardinality of N equals, c, the cardinality 

of continuum. Let g : N — > (—1,0) denote a natural 1-1 correspondence. Now, there exists 

an injection fi : P{R) -> P, defined by f\(A) = g(A), when A e N; f\(A) = K, when eq(|) 

admits infinitely differentiable solutions in A; and fi(A) = K p , otherwise. Conversely, one also 

finds an injection f 2 : P —> P(R), defined by f2(K) = {K}, but f2(K p ) = {K, B}, where 

3 From now on, we denote by <f> this stochastic time-like feature of the golden mean. By the symbol v we, 
however, continue to refer to the usual 'non-random' number ^T 1 ■ 
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B = (0, 1) is the interval where K p is defined. Hence, by the Schroeder-Bernstien theorem [7], 
the cardinality of the physical set P equals 2 C . Note that the mapping g renders all possible 
discrete subsets of P unphysical, transferring them to (-1,0). Further, eq(|8]) admits infinitely 
differentiate solutions, over a given, apriori scale denoted as 77, provided the possibility of 
random scale dependence via inversions is neglected. Consequently, the set R is realized in P 
in an approximate sense. Moreover, there exist an uncountably more nontrivial elements in the 
physical set in comparison to the real set. This means, in particular, that there exist physical 
numbers ( in the form of fluctuations) which are infinitely smaller than any nonzero ordinary 
real numbers. By inversion, one then concludes that there exist infinitely large physical numbers 
greater than any real number. Clearly, all the set theoretic results, valid in R , get carried over 
to the physical set. In particular, the physical set is partially ordered, when k p < K p means 
k < K. 

Now to construct a nontrivial infinitely small fluctuation, let us recall that a physical number 
is associated with an accuracy limit. Let k p = rk(r]~ l )r) be determined with an accuracy 10~ m . 
Let sup[l/fc] ~ iQ( m + s ) , Here, [.] denotes the greatest integer function, and the supremum is 
defined on a bounded interval of 77. Then, as m — > 00, s large, but finite, k p = 0, in the ordinary 
sense, for finite 77. But in the sense of a 'physical limit' 771 = 10~ m i] — > 1, m — > 00, 77 — > 00, one 
gets an arbitrarily small random number, by allowing s to assume larger and larger values, but 
never exactly allowing s = 00. Because of randomness, k p is not an ordinary real number. The 
physical limit tells that the ideal condition of infinite accuracy is realised only in an infinitely 
distant time. In other words, any measurement process over a finite period of time can achieve 
only a finite degree of accuracy. We note that the physical set P has a structure analogous to 
the nonstandard real number set [5]. A more detail investigation of the relationship between 
the two approaches will be considered separately. 

Next we show how such an infinitely small fluctuation k p can have nontrivial influence at 
the level of the ordinary scale, thereby re-deriving, in an alternative way, the small scale time 
evolution in the golden mean v. Let us rewrite the physical time t in the form t = 77^(77) , <p(i]) = 
1 + kpcj)^ 1 ), rji = krj. Since t must be an exact solution of eq(||), one obtains k P d ^ = 0. 
It thus follows that in the ordinary scale of r] : = 0, and hence <\> is an ordinary constant 
(having no evolution). However, this conclusion may not necessarily be true at the level of 
the smaller scale 771. In fact, as i]\ grows to order 0(1), d ^ , need not vanish, since k p , 
though nonzero, becomes vanishingly small provided d ^ ~ O(l)- To verify this, note that 
0(7/) — 1 = 77i(0(7/]~ 1 )) ((p < 1 as long as r\\ < 1), 771 ~ 0, so that in the limit of 771 = 1 + 771 — ► 1, 
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^ = ^) (9) 
dm 771 

which is of order 0(1), and one retrieves a self similar replica of eq(H), over the scale 771. Note 
that, as 771 — > 1~, <p behaves as a small scale linear time so that <fi— 1 « — d(p. But as rji- — > 771+, 
by inversion, (-) sign cancels so as to reproduce eq(P). It is now easy to verify the continuity 
of the first derivative at rji = 1, since ~ ~j^~=r = dtn?)i = ^he continuity of the second 
derivative can not, however, be maintained because of the sign difference injected by inversion. 

Note that the above analysis leads us to the fact, by yet another route, that the golden mean 
v = (f) — 1, must have small scale intrinsic time evolution. To state more clearly how a time 
sense gets attached to v, let us note that the intrinsic time must flow, in the neighbourhood 
of T)(— 1), in the form <j){rj) = 1 + k p (j){j]^ 1 ), till a random inversion near r\ = k~ x carries it to 
the smaller scale 771, to assume the form 0(77) = 1 + /c(0(?7i)) _1 . As a consequence, a (linear) 
time sense, self similar to the scale rj is generated over the scale 771 , which persists upto order 
f]i ~ k^ 1 , preparing it for yet another replication on the second generation scale 772 and so on, 
leading to an infinite continued fraction = 1 + [k; k;k; . . .] representation of the intrinsic time 
flow, as the zeroth generation ordinary time rj — > 00, exploring longer and longer scales. Letting 
k — 1, one gets the golden mean flow of time. Note that the intrinsic sense of time is derived 
from the infinitely slow random unfolding of cascaded scales hidden in the form of an infinite 
continued fraction. Among all these time-like continued fractions <pk, the golden mean continued 
fraction is distinguished by its slowest possible convergence (unfolding) rate. We now show that 
<pk ~ 1^0as7/^oo. Note that after nth inversion, 77 reaches the scale r] n . Let fi M denote 
the smallest infinitely large natural number, exceeding all real numbers, in the physical set P, 
so that r/oo = k n °°r]. One exhausts all possible ordinary scales by utilising all the available 
number of inversions so that the final replication leads to 0(t7oo) = 1 + ^ \ , yielding v. It thus 
follows that although, in general, time may flow following the steps of fa, this scale dependent 
flow can in fact continue at most upto a finitely many scales (in the context of physical time). 
The intrinsic flow, exploring smaller and smaller scales at slower and slower rates, as the 
ordinary zeroth generation scale 77 — > 00, would finally converge to the slowest possible golden 
mean flow. Note that once the golden mean flow is reached, no further scale replication is 
allowed, since all the higher order scales would be physically indistinguishable (7700+1 — Voo), 
from the stand point of the ordinary scale 77 ~ 1. However, because of self similarity, any scale 
77„ could be considered as the zeroth generation scale, the process of exploring the golden mean 
flow of time would remain unaltered. 

We close this section with a few more remarks, highlighting a number of important features 
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of the golden mean. 

1. The golden mean time sense is intrinsic, since it is independent of an apriori time. Note 
that 0(0(.)) = 0(.), where (.) indicates that could be a function of any t n . In fact, the 
equation means = 0(0), thus eliminating the Newtonian external time rj. (Recall that to 
avoid any confusion, we choose to indicate the slow scale dependent (logarithmic) variation 
in the golden mean, while v~ x denotes the usual constant value of it, over a well defined given 
scale.) However, an approximate, 'coarse-grained' Newtonian sense of time is realised over 
a scale, when the small scale fluctuations, because of local inversions, are ignored. In the 
following, we denote by 0(0) the set (equivalence class) of all possible intrinsic variables of the 
form {(f)(t n )}. 

2. So far, we have not spelled out the form of the random variable r. A probability 
distribution satisfying the constraint < r n , < r >= r can be written as follows. Let the 

sample space of r be {ro +1 }, s = 0, 1, 2, The corresponding discrete probability function is 

defined by p(r) = e~ r o 1 ^,-, m being a positive integer. The distribution is Poisson-like, but not 
exactly the same, because of the special sample space, generating the infinite set of scales. Let 
r < 1, so that r™ « 1, for a suitably large m. Thus < r n >= rQp(ro),p(r ) = e^^o -1 ) m 1. 
Further, any number r can be written as r = u a , a = log^ro- Thus all the free parameters in 
the definition of the fractal time t is determined self consistently in terms of v. The parameter 
k — 1, since the scale factor is included in the sample space. One thus fixes A in eq(l) and (6) 
as A = r . 

3. The golden mean v remains constant over a scale 77, but, nevertheless, enjoys intrinsic 
randomness. Because of the equality v — (1 — z/)£z/ 2n , v is realized as the expectation value 
of a r.v. r with sample space {v n } of scales, the corresponding probabilities being {p(r)} = 
{(1 — v)v n }. Note that v here denotes the exact value represented as Let us consider 
the sample space {S m+n } of yet another r.v. r„, for a sufficiently large m, where S n is the 
nth approximants of the above series. The corresponding probability distribution is p{r v ). 
Consequently, r v realizes higher precision values of v with lower and lower probabilities. Clearly, 
v = r v . This proof, being purely of number theoretic origin, tells that measurements in a 
physical universe are inherently uncertain, since any measurable quantity is a multiple of v. 
Note that the evaluation of the exact v needs a measurement process defined over an infinite 
amount of time, thus making it physically impossible. The choice of the probability distributions 
may, however, appear arbitrary. In Sec. 6. 2, we show how these discrete distributions relate 
naturally to a universal density function for the class of fluctuations considered in this work. 
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5 Applications 



In the following applications of the fractal time, we show how the formalism yields (i) a natural 
resolution to the generic 1/f signal problem [1], and (ii) the universal probability density 
observed in the Bramwell- Holdsworth- Pinton (BHP) fluctuations [2]. 

To begin with let us first reconsider eq(|lj) in the light of fractal time t. Recalling t = 
77(1 + r0(r/ -1 )), dt — (1 + r0)dr/, since ^ = 0. The equation (f7|) in the 'physical' time t is 
re-expressed in the Newtonian time as 



n t 

— = (l + r0(r / ~ 1 ))x (10) 
d?7 

An exact solution of this equation, eq(^), as an implicit function of 77 has already been obtained. 
An analogous form of the solution is now written as x = e^ 01 ^ a{rf) = r J r]d(j){r]~ l ) , when 
eq([|) is used. Note that the infinite series of scales T n , with alternating signs, provides an 
explicit representation of the intrinsic time variable 0, when the system is allowed to evolve over 
all the available scales in an infinite period of the ordinary time 77. The exponent a thus leads 
to the exponent fi of eq(|jD, when the integration is performed successively over scales between 
two consecutive moments of inversion. It thus follows that the method of Sec. 3 actually yields 
the fractal solutions to eq ([?]), when k is identified with one of the scales v n . Moreover, a 
solution of eq(|]) when k is treated as a system variable can be obtained as x = e 7?+ ^ KJ?2_/3 ^ , 
where f3{rj) = r j (j] + 2krj 2 )d(f)(r]~ 1 )) , neglecting 0(r 2 ) corrections. Note that both the solutions 
indicate almost identical late time power law fluctuations. 

5.1 1/f spectrum 

To calculate the spectrum of the stochastic fluctuations, present universally in linear equations 
of the form eq(|]), we need to estimate the late time asymptotic form of the corresponding two- 
point correlation function C(r]) =< x(0)xf(r]) >=< Xf(rj) >, since x(0) = 1. Assuming that 
|A| ~ v n (say), it follows that C(rj) ~ 77™^, and hence the spectrum has the form S(f) ~ 1// 1=FM , 
\i being the slowly varying function of Sec. 3. Clearly this should be the generic form of the 
spectrum for a general dynamical system of the form 

n t 

— = h(t)x (11) 

where the time dependence in h may have nonlinear influences: h(t) = h(t,x). However, the 
nature of explicit nonlinearity in a system is expected to get reflected in the exponent /1. The 
generic logarithmic correction in \i provides significant insights into the late time features of 
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the dynamics, which might get revealed in a time series over a number of different scales. One 
example is treated below. 

5.2 Universal probability function 

Recently, a universal pattern of self similar fluctuations have been reported to occur in many 
Natural processes. Subsequently, a generic probability density function (PDF) of the form 

Pit) = Ke at - aet (12) 

is shown to be respected by the underlying dynamics [2] , of apparently unrelated systems. Here, 
t is a relevant fluctuating variable. Let us note that the solution eq(|]) denotes the universal 
fluctuation pattern, at least, for those Natural processes which satisfy an equation of the form 
eq(|ll~D. We now show how the above PDF is naturally realized for this universal fluctuation. 

Let us begin by noting that the infinite alternating series of scale dependent terms in eq(§) 
gives the complete fluctuation spectrum Xf of a linear system, which is self similar over all 
these scales and is revealed over an infinite period of the ordinary time i]. However, the first 
two terms in the series are sufficient to capture the generic features of the fluctuation Xf , 
because the scale generating r.v. r (c.f., remark 2, Sec. 4) induces a higher order (stochastic) 
scale dependence on each of the scales T n , thereby inscribing a complete replica, of the total 
fluctuation Xf over the scale T\ (say), even in a finite period of rj. Further, the (-) sign between 
two consecutive scales T% and T 2 is a nontrivial signature of inversion. Thus it suffices for us 
to consider a renormalised fluctuation of the form Xf oc e _r< - Tl_T2 - ) (A = ro), for a finite 77. 
Accordingly, the generic PDF corresponding to Xf should be identical with the same for Xf. 

Let, for definiteness, that the random scales r = {z/ n-1 } be distributed with probabilities 
p(r) introduced in remark 3, Sec. 4. Because of the logarithmic scale dependence T 2 = ln(l + 
u(Ti — 1)), one gets Xf oc e r( - T2 ~ ( - 1+u ' ,eT2 \ Dropping the higher order scale dependence introduced 
by the v dependent factor, we get finally Xf oc e r ( T2 ~ eT2 \ Now in a practical situation, the 
fluctuation Xf is represented in the form of a time series record over a period of time. The 
probability that a randomly sampled observation Xf n is drawn from the scale r n , in the sample 
space of r, now equals the product of the probability of selecting the scale r n and the conditional 
probability that the observation actually comes from the said scale, given that the scale has 
been chosen already. But the conditional probability is nothing but the correlation function 
C(T 2 ) =< x fn (T 2 )x f (0) >=< Xf n (T 2 ) >, 2/(0) = 1, x fn oc e r "( T2 " eT2 ). Thus the probability of 
drawing a random sample from the scale r n equals P{r n ) oc p(r n )e rn( - T2 ~ eT2 \ Hence the grand 
universal probability that the time series reveals the whole spectrum of fluctuations over all 
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possible scales is obtained as 

P u oc ££°P(r n ) = (1 - u)T^u n e {l+v)vn{T2 - eT2) (13) 

Clearly, apart from a multiplicative factor, which could be fixed from the normalisation condi- 
tion, the zeroth order term of this infinite series representation of the universal PDF agrees well 
with that in ref.[2], which was obtained from an approximate argument (the factor a = 7r/2 
in the exponential gets replaced here by 1 + v). The higher order terms in the infinite series 
represent corrections, which are likely to give more accurate fits of the time series records for 
Natural processes, as noted in ref.[2]. 

Now to explain the reason of the matching, we note that the factor 1 + v in eq([T3|) actually 
realises T2 as T\(=. (1 + v)T<i) and T\ as rj. In the case of an explicitly nonlinear system given 
by eq(|TTD with nonlinearity coupling k ~ 1, the system experiences fluctuations at a time 
rj ~ 1. The scale T\ then corresponds to the 2nd generation fluctuation in a corresponding time 
series record. It thus follows that the zeroth order PDF would have the form e Tl ~ eTl . Let us 
note that the moments of the model fluctuating variate in ref . [2] possesses the generic property 
< r n >oc Tq. However, the corresponding PDF is obtained from a quantum field theoretic 
consideration of a critical magnetic model which is based purely in the Newtonian time frame. 
Now, the relationship between an 'ordinary' dynamical variable Qo, following an evolutionary 
equation of the form eq([ll|), but in the ordinary time, and the corresponding 'physical' variable 
Q p is given by Qo oc Qp 1 , since in the logarithmic scale \nt = (1 + v) In?]. An application of this 
conversion rule thus leads to 

P BHP oc e ( 1+ ^" eTl ) (14) 

which completes the derivation of the BHP probability function. 

To proceed further, let us now re-derive the generic PDF from an alternative method. This 
will reveal a subtle relationship between the universal PDF and the set of discrete distributions 
we have chosen as examples. Let r] n = u n i], where r] e (0, 00) be nonrandom, and v denote an 
approximate value of the golden mean. The scales i] n are random, because of the randomness in 
v and is assumed to follow the distribution of remark 2, Sec. 4 (m = 1), so that 771 is realized with 
probability Pi = e~ v v. Since v is approximate, it will now undergo evolution in the physical 
set following, for instance, the scale free representation v v = ucj)^ 1 ) = vj{\ + </>(^ii )) = 
f/(l + /{I + ^(Vii)})' e ^ c - We may assume that the linear sense of time generated by <\> at 
the scale of i] is denoted by 77 itself. As 77 flows from 77 w slowly, the intrinsic evolution in 
v splits the scale 771 into an infinite set of tiny scales rji n = v n rji, each of which will further 
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undergo finer levels of subdivisions, and so on. The stochastic evolution of v thus fractures the 
scale r\i at the neighbourhood of every point, thereby raising it to the level of a continuous r.v. 
with the PDF Fifa) oc e~ Hi ^ ) 4>(r] 1 1 ), a Gamma distribution. We note that El Naschie [8] 
has already used a Gamma distribution to derive the Hausdorff dimension of the fluctuating 
(Cantorian) spacetime as 4 + u 3 . He has also indicated how the mass spectrum of all known 
elementary particles could be determined using the golden mean v. (The uncertainty exponent 
of the fractal time is obtained independently by the author as v [10].) 

All the higher order scales will similarly undergo infinitesimal fluctuations and hence finally 
be distributed following the above Gamma density function, since = v n 4>{rj~ l ) . Note that all 
these (f){f]n l ) functions correspond to different scale dependent realizations of the same universal 
function 0(0) in the limit of infinite time. Letting z/0 = e^, we reproduce the universal PDF 
for scale free (self similar) fluctuations, P(<j>) oc e* -60 . But is again a realisation of 0(0), by 
eq(||) and remark 1, Sec. 4. 

6 Concluding remarks 

The extension of the real set to the physical set provides a dynamical representation of the 
number system, each member of which is associated with an equivalence class of fluctuating 
elements separated by infinitesimally small scales. The fundamentally stochastic nature of the 
golden mean number renders accurate measurements in the physical set impossible. Conse- 
quently, the physical universe based on the physical set would consists of intrinsic changes and 
fluctuations. It is hard to imagine any fundamentally constant physical quantity in this uni- 
verse. The real number set is an incomplete realization of the physical set, when the possibility 
of infinitesimal changes by local inversions are ignored. However, because of the new exact 
class of solutions, P = R. In the midst of all these changes and approximations, there exists, 
however, one symbol of perfection in the form of the golden mean equation 0(0) 2 + 0(0) = 1, 
being engraved fundamentally in the formalism of the SL(2,R) Calculus. A more detailed, in 
depth analysis of this Calculus, will be presented in a subsequent paper, where the status of 
well known theorems such as the Picard's existence and uniqueness theorem will be examined. 
Let us only remark here that the present class of solutions are not in contradiction with the 
Picard's theorem, the scope of which gets extended in the SL(2,R) formalism. The present 
formalism is likely to initiate a new approach in understanding the origin and dynamics of the 
nonlinear phenomena in Nature. 
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